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ABSTRACT 



Context. Hot- Jupiter planets must form at large separations from their host stars where the temperatures are cool enough for their 
cores to condense. They then migrate inwards to their current observed orbital separations. Different theories of how this migration 
occurs lead to varying distributions of orbital eccentricity and the alignment between the rotation axis of the star and the orbital axis 
of the planet. 

Aims. The spin-orbit alignment of a transiting system is revealed via the Rossiter-McLaughlin effect, which is the anomaly present in 
the radial velocity measurements of the rotating star during transit due to the planet blocking some of the starlight. In this paper we 
aim to measure the spin-orbit alignment of the WASP-3 system via a new way of analysing the Rossiter-McLaughlin observations. 
Methods. We apply a new tomographic method for analysing the time variable asymmetry of stellar line profiles caused by the 
Rossiter-McLaughlin effect. This new method eliminates the systematic error inherent in previous methods used to analyse the effect. 
Results. We find a value for the projected stellar spin rate of vsin; = 13.9 + 0.03 km s -1 which is in agreement with previous 
measurements but has a much higher precision. The system is found to be well aligned, with A = 5*f° which favours an evolutionary 
history for WASP-3b involving migration through tidal interactions with a protoplanetary disc. From comparison with isochrones we 
put an upper limit on the age of the star of 2 Gyr. 

Key words, planetary systems - eclipses - Techniques: spectroscopic - Line: profiles 



1. Introduction 

Since the confirmation in 1995 of the first extrasolar planet 
orbiting a main sequence star dMavor & Ouelozl 1 1995b it has 
become evident that gas-giant planets with orbital separations 
smaller than that of Mercury are quite common. In fact, these 
Hot Jupiters account for around 20% of the approximately 45CQ 
exoplanets discovered to date, though their current prevalence 
in transit searches is mostly due to observational bias and re- 
cent results from the Kepler Mis sion suggest that the majority 
of ex oplanets are of lower mass (Borucki & for the Kepler Team! 
2010). It has been shown that these planets must form farther 
out, b eyond what is known as the 'snow line' dSasselov & Lecad 
2000 ), and the n migrate inwards to their current locations 
dLin et ail 1 19961) . If this is the case then through what process 
does this migration occur and how can we gain information on 
this from current observations? 

There are three main competing theories of Hot-Jupiter mi- 
gration. The first states that migration occurs due to tidal in- 
teractions between the planet and the protoplanetary disc caus- 
ing some of the planet's angular momentum to be lo st to 
the disc (Gold reich & TremainellT980l: iNelson et alj|2000 !). The 



For up-to-date exoplanet statistics see http://www.exoplanet.eu 



second proposes that gravitational scattering between multi- 
ple planets in a system can cause one planet to migrate in- 
wards at the expense of the other bei n g shot out of the sys- 
tem dWeidenschilling & Marzaril [19961: Ichatteriee etail 120081: 
lJuric & Tremainell2008l) .~ 

These separate scenarios should lead to quite different post- 
migration states. Tidal interactions between the planet and pro- 
toplanetary disc will result in a system where the orbital axis 
of the planet is well aligned with the rotation axis of the host 
star. It would also lead to highly circularised orbits. The planet- 
planet scattering process on the other hand would result in plan- 
ets which have l arger orbital eccentricities and are not necessar- 
ily well aligned ( Weidenschilling & Marzari 1996). 

The third proposed process involved in Hot Jupiter migra- 
tion is the Kozai mechanism. It involves oscillations between 
the spin-orbit alignment angle and the orbital eccentricity of 
the planet due to the presence of another (outer) planet in the 
system or a binary companion of the host star dKozai| [l962; 
IWu & Murravll2003UFabrvckv & Tremainell2007l) . During these 
oscillations the value - e 2 cos i is conserved, where e is the 
planet's orbital eccentricity a nd ; is its inclination to the orbital 
plane of the other objects. Nagas awa et al.l d2008l) suggested a 
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combination of all three processes could be responsible for the 
migration of Hot Jupiters. 

Measuring the alignments of the known Hot-Jupiter systems 
helps shed some light on the migration process and refine mod- 
els of system evolution. The misalignment angle A between the 
planet's orbital axis and the rotation axis of the host star can be 
determined via the Rossiter-McLaughlin effect dRossiterill924t 
McLaug hHrIT924h . This is the radial velocity anomaly observed 
during transit due to the planet blocking some of the starlight. 
A planet in a prograde orbit would first block some blue-shifted 
light, from the half of the star which is rotating towards us, caus- 
ing an anomalous red-shift in the star's radial velocity. Then once 
it has crossed the stellar rotation axis it will block light from the 
receding half of the disc causing an anomalous blue-shift. The 
form of this radial velocity anomaly allows us to measure the 
projec ted stellar rotation ra te and the spin-orbit misalignment 
angle (iGaudi & Winnll2007l) . 

lOhta et alJ d2005l) and lGimenezl (2006) derived detailed ex- 
pressions for measuring the anomalous shift in the line cen- 
troid during transit. Using the radial velocity shift to deter- 
mine the parameters related to the Rossiter-McLaughlin effect 
can be a problem when the spectrograph can resolve the stellar 
line profile. The data pipelines for instruments such as SOPHIE 
and HARPS calculate RVs by fitting a Gaussian profile to the 
cross-correlation function (CCF) in order to measure the shift 
in line centroid. However, as the planet is blocking some of the 
light from the star this shows up on the line profile as a trav- 
elling 'bump' of width equal to that of the local non-rotating 
line profile. This means that the CCF is now asymmetric and 
time-varying during transit, introducin g a systematic dev iation 
to measurements of the radial velocity (I Winn et al J 120 05). This 
error becomes greater for more rapidly rotating host stars as 
their line profiles exhibit higher degrees of rotational broaden- 
ing dHirano et al . 200§. 

In more recent studies IWinn et all d2005l I2006L 120071) tried 
to account for this by building models of the out of transit 
line profile and the light blocked by the planet. They included 
these into their analysis of the line-spread function and pro- 
duced semi-empirical corrections to the model radial velocities. 
However, their method does not account for the problem entirely 
as their analysis of HD 189733b shows a clear patte rn of corre- 
lated residuals in the radial velocities during transit dWinn et al.l 
2006). A similar pattern in the radia l veloc ity residuals of HD 
189733b was f ound by iTriaud et alJ d2009l) who discuss their 
cause in detail. ISimpson et alJ d2010h also found a similar cor- 
relation pattern when analys i ng the WASP-3 system using the 
expressi ons fro m Ohta et al. (2005) with the corrections devel- 
oped bv lHirano et alJ (|2009l). In order to eliminate t he need for 
empirical corrections ICollier Ca meron et alJ d200 9a) developed 
a new method that involves decomposing the CCF profile into its 
various components, namely the limb-darkened rotation profile, 
gaussian average line profile and the travelling signature caused 
by the transiting planet. 

To gain information on the distribution of spin-orbit mis- 
alignment angles, the WASP (Wide Angle Search for Planets) 
consortium has been performing follow-up spectrographic 
observations of the transits of the known WASP planets. 
Observations have been made using the HARPS spectrograph 
on the ESO 3.6m telescope at La Silla and the SOPHIE spec- 
trograph on the 1.93m telescope at the Observatoire de Haute- 
Provence. WASP-3b is the third planet d iscovered through the 
SuperWASP project dPoUacco et al.ll2008l) . It is a Hot Jupiter or- 
biting a main sequence star of spectral type F7-8V. 



In this study we implement the method of 
ICollier Cameron et al.l d2009al) to remove the inherent er- 
ror present in the Gaussian-fitting method by instead focusing 
on the light removed by the planet. We model all the components 
of the line-spread function and fit to spectral observations of 
WASP-3b in order to track the trajectory of the missing light 
component as it crosses the line profile during transit. A model 
of the out-of-transit profile is subtracted from the data leaving 
only the signature of the blocked light. The trajectory of this 
feature is used to derive values for the projected stellar rotation 
rate, the impact parameter and the spin-orbit misalignment. 

This method does not rely on high precision radial velocity 
measurements, making it uniquely usef ul for detecting planets 
orbiti ng rapid-rotating early-type stars dCollier Cameron et al.l 
2010). Many of these stars which have reliable photometric data 
have been dropped from radial velocity follow-up observations 
due to their rotation rate and line-poor spectra. 

2. Observations and analysis 

We reanalyse data presented in Simps on et al. ( 20101) w ho used 
the SOPHIE echelle spectrograph (Bouchv et alJ 20091) on the 
1 .93m telescope at the Observatoire de Haute-Provence to take 
26 observations before, during and after the transit of WASP- 
3b on the night of September 30th 2008. A major advantage 
of this new method is that it does not require high precision 
radial velocity measurements, therefore the spectrograph was 
used in high efficiency mode (R = 40000). The exposure times 
for the observations ranged from 300 - 1800 seconds to en- 
sure a constant signal-to-noise ration of 35. In total 137 minutes 
of observations were taken during transit and 130 minutes out 
of transit. CCFs were computed using the automated SOPHIE 
data-reduction pipeline which is adapted from the HARPS data- 
reduction software. WASP-3 is of spectral type F7-8V, so the 
weighted mask function used for the cross-correlation was for a 
G2V star as this was the closest type offered by the pipeline. For 
a more detailed ex planation of the obse rvations and data reduc- 
tion procedure see Sim pson et al.l d2010t) . 



2.1. Modelling the observations 

The following sections outline the key points of the new method 
used to analyse th e CCF data. The modelling p rocess is de- 
scribed in detail bv lCollier Cameron et al.l (l2009al) . 

First we need to build a model of the out-of-transit CCF. 
This is done by convolving a Gaussian representing the the local 
line profile at any point on the surface of the star with a limb- 
darkened rotation profile. For the local line profile we use 



g(x) 



1 

e & 



For the limb-darkened rotation profile we use the equation 



/(*) 



6((1 -m)V1 -x 2 
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n(u - 3) 



(1) 



(2) 



assuming a linear limb-darkening model where u is the limb- 
darkening coefficient and — 1 < x < 1. We adopted a value of 
u = 0.69 taken from the Claret (2004) tables for the g' filter. 
This value corresponds to a star with T e ff = 6500 K, logg* = 
4.5 and [M/H] = which best describe the values for WASP-3 
calculated through analysis of the SOPHIE spectra and presented 
in the discovery paper (T eff = 6400± 100K, log g, = 4.5+0.05). 



2 



G. R. M. Miller et al.: The Doppler Shadow of WASP-3b 



The convolution of these two functions is given by 



2.3. Markov chain Monte Carlo technique 



h(x) 



1 

I 



f(z)g(x - z)dz 



(3) 



which is calculated by numerical integration. We need to 
shift the model CCF to account for the fact that the CCFs are 
computed in the velocity frame of the solar system barycentre. 
To do this we compute 



Xjj = Vij - (K(e cos a) + cos(vy + lo)) + y) 



(4) 



where v, y is the velocity of pixel i in the barycentric frame, 
Vj is the true anomaly at the time of the y'th observation, w is 
the argument of periastron, e is the orbital eccentricity, K is the 
radial velocity amplitude and y is the systemic centre-of-mass 
velocity. 

At any moment in time the position of the planet on the plane 
of the sky is given by the co-ordinates 



x p — r sin(v + co — n/2) 
z„-r cos(v + oj - n/2) cos i 



(5) 
(6) 



where r is the instantaneous distance of the planet from the 
star and i is the inclination of the orbital axis to the line-of-sight. 
The perpendicular distance of the planet from the stellar rotation 
axis in units of R, is now 



u p — x p cos A — z p sin A 



where A - 



Voi-bit 



(7) 

and <p is the position angle in the 



plane of the sky dWinn et al.ll2005l) . Combining the model out- 
of-transit CCF with the model of the missing starlight gives us 



Pij = h(xjj) +Pg(xij - Up) 



(8) 



Here the term A(jCy) is the model of the out-of-transit stellar 
CCF. The term /3g(xij - u p ) represents the travelling Gaussian 
signature caused by the missing starlight, where ft is the frac- 
tion of starlight blocked by the planet during the total part of the 
eclipse. Fig. 2 shows the original CCFs and the resulting residu- 
als when first the model CCFs h(xjj) are subtracted off leaving 
the signature of the missing starlight and secondly the overall 
residuals when the complete model py is subtracted. 

2.2. Fitting the model 

In order to orthogonalise the data and the model we subtract their 
optimal mean values using inverse-variance weights w,j = 1 [erf-: 
The goodness of fit of the model to the observed data is then 
calculated using the^- 2 statistic 



X 



(Oh 



(9) 



where a, is the optimal average of the residual spectra and A 
is a multiplicative constant calculated via optimal scaling. 

The spectral resolving power of the SOPHIE spectrograph 
is R=75000 which gives a velocity resolution of 4 km s . 
However, the CCFs are produced with velocity increments of 
0.5 km s _I meaning that the errors on the data from neighbour- 
ing pixels are correlated. In order to make sure our data points 
were statistically independent we binned them by a factor of 8 
before computing^ 2 . 



The model parameters were calculated using a Markov chain 
Monte Carlo (MCMC) method. The MCMC code used for this 
study is a hybrid of the code previously used to calculate the 
parameters of the WASP systems from photo metric and spectro- 
scopic data s ets dCollier Cameron et al.| 20071) . and a new code 
developed bv lCollier Cameron et al. (2009a) specifically for RM 
analysis. Therefore at the same time as calculating the parame- 
ters from the RM effect we recalculated all the MCMC fitting pa- 
rameters for the system. A compar ison of our results wit h those 
of the initial MCMC analysis from lPollacco et al.l d2008l) can be 
seen in Table 1 . The co de includes the mass and radius calibra- 
tion dTorres et al.ll2010|) which was rece ntly implemented in the 
M CMC analysis bv|Enoch et alj d2010t) . The method described 
bv lTorres et all (12010) " is used to derive the stellar mass and ra- 
di us from polynomial functions of T e ff, logg* and [Fe/H]. As 
in ICollier Cameron et al.l (l2009al) we replaced the width of the 
local non-rotating profile s with vqcf, the FWHM of the CCF in 
an attempt to avoid correlated pairs of parameters. 



v C cf = ^(vsinO 2 + v 2 



(10) 



where v g is the FWHM of the Gaussian representing the local 
stellar and instrumental line profile. 



v g = 25V sin i 



(ID 



At each step in the chain the current values of the four pa- 
rameters are altered by a Gaussian perturbation and then the 
goodness-of-fit is recalculated. For example, the proposed next 
step for A is 



h = h-x + fcr A G(0, 1) 



(12) 



where / is a scale factor of order unity and G(0, 1) is a ran- 
dom number drawn from a Gaussian distribution of zero mean 
and unit variance. Steps are accepted or rejected in accordance 
with the Metropolis-Hastings algorithm. After each proposed 
step A^ 2 = x\ ~ x\-\ I s calculated. If A;^ 2 < then the pro- 
posed step is accepted. If A% 2 > then the step is accepted with 
a probability of e~ Ax ^ 2 . Each successful step is recorded to the 
MCMC output. If a step is rejected then the previous accepted 
step is recorded again. Varying the scale factor / changes the 
acceptance rate of the MCMC. We ran the MCMC with an ini- 
tial burn-in phase of 1000 steps. After the first burn-in period we 
re-evaluated our estimates of the variances on the binned CCF 
data. We then ran the chain for another 100 steps in order to 
re-evaluate the variances on the four fitting parameters from the 
chains. This was followed by a final production run of 10,000 
steps where / was fixed at a value of 0.5 as this was found to 
return the desired acceptance rate of 25%. 

2.4. Photometric and spectroscopic datasets 

In total we used 8 photometric and 2 radial velocity datasets 
in the MCMC analysis. In addition to the photometric datasets 
analysed in the discovery paper we included 2 additional sets of 
photometry taken with the RISE instrument on the 2m Liverpool 
Telesc ope at Observatorio del Roque de Los Muchachos, La 
Palma (Gibson et al. 2008). The radial velocity data comprised 
the 26 observations described earlier to target the RM effect, 
and the 6 out-of-transit observations taken in July and August 
2007 also using the SOPHIE spectrograph as described by 
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Fig. 1. Upper panel: Phase-folded plot of the 6 out-of-transit radial velocity measurements which are not affected by any time- 
varying asymmetry of the line-profile. The out-of-transit RV fit was calculated using the velocity semi-amplitude, orbital eccen- 
tricity, arguement of periastron and true anomaly. The position of the planet over the stellar disc, ratio of star/planet radii, impact 
parameter and non-linear limb darkening coefficients were used to model the RM anomaly during transit. Lower panel: Phase-folded 
plot of all 8 sets of photometric data analysed in this study. 



Polla cco et al.l d2008l) . Figure 3 shows phase-folded plots of the 
photometric data and the out-of-transit radial velocity measure- 
ments. 



3. Results 

The resulting values for the system parameters can been seen 
in Table 1 alongside those taken from the discovery paper, both 
with their l<x errors. With the exception of the impact param- 
eter, all values are in agreement within their l<x uncertainty 
ranges. Throu gh examinat i on of its relation with the stellar sur- 
face gravitv lPollacco et al.l d2008l) suggest that the impact param- 
eter should lie between 0.4 and 0.6 which is in good agreement 
with the value found in this study of b — 0.38^ ^. Using our 
method the impact parameter can be more closely constrained, 
as the properties of the streak of 'missing light' also give us a 
measure of the latitudes of ingress and egress. 

One major success of this study is that sensible values in 
agreement with previous work are produce d without having to 
fix any of the system parameters in place. IPollacco et al] (f2008) 
showed that in order to reconcile the MCMC analysis with spec- 
troscopic diag nostics log ff, must be between 4.25 and 4.35. 
Recently Eno ch et alj ((2010) attempted to reproduce the results 
from the WASP-3 discovery paper incorporating the Torres mass 
calibration into the MCMC analysis and found they also needed 
to fix log g» to the value determined by spectroscopic analysis. In 
this study we left log g* as a floating variable and found its value 



to be logg, = 4.33^q q3 which lies within the range suggested in 
the discovery paper. 

The results produced for the RM parameters are presented in 
Table 2 along with their one sigma errors. The value we obtained 
for the projected stellar rotation rate v sin i = 13.9 + 0.03 km s _1 
is in good agreement with the value of v sin i — 13.4+ 1 .5 km s _1 
derived by the analysis of the SOPHIE s pectroscopy a s presented 
in the WASP-3 discovery paper (IPollacco et al.l l2008). However, 
our result shows a much greater level of precision. This is be- 
cause uncertainties on previous estimates of v sin i are removed 
by the fact that we can measure the FWHM of the intrinsic pro- 
file directly. Our analysis finds the projected spin-orbit misalign- 
ment angle to be A = 5*5°, which is almost indistinguishable 
from zero. 

ISimpson et al] (1201 Oh recently performed a Rossiter- 
McLaughlin effect analysis of the WASP-3 system with the same 
data u sed in this study. They used the method of iHirano et al.l 
(2009) to account for the aforementioned systematic error en- 
countered when trying to calculate vsin/ by applying the 
lOhta et al.l d2005l) method. Their value for v sin; of 15. 7+ J \ km 
s _1 is similar to our value but still slightly larger. For A they de- 
rived a value of 13 + ^° whic h agrees with our result. In addition 
to this Tripathi 'et al.l (1201 Ol) performed an analysis of a separate 
set of data and retrieved values that are in good agreement with 
those found in this study, as can be seen in Table 2. 

Using our new, more accurate values for the system pa- 
rameters we were able to plot the position of WASP-3 on the 
R/M l/3 vs T eff plane (see Fig. 3). From the MCMC analy- 
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Table 1. Results of MCMC analysis compared to those presented in the literatures. 



Parameter 


Symbol 


Pollacco et al. (2007) 


Gibson et al. (2008) 


This study 


Units 


Transit epoch (JD-2450000) 
Orbital period 
Planet/star radius ratio 
Transit duration 
Impact parameter 


T Q 
P 

(Rp/R*) 2 
t T 
b 


4143.8503+S 

1 S4.fiR^4. +(,on(l1 " 12 

0106 +,io002 
1H0+™ 3 
505 +0 051 

u,JUJ -0 166 


4605 5592+ 00002 

1 S4.fiR^^ +0 ' (ll),l0,c 
i.otuojj_ 0000002 

0103+™ 002 
1 1 47+O.0TO 
u.nt/ _gnnp5 

0.448!™ 


3992 4101+ 0004 

1 C/l AOQ7+0.000001 
i.S4t38J/ 0Q0000] 

0105 +0 0002 
0-1126™ 
0.38_*J 


days 
days 

days 
R. 


Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semi-major axis 
Orbital inclination 


K s 

r 

a 
i 


9 S 1 9+0 0079 

U,ZJ1 -00108 

S /LSR7+000I3 
-3.4S8/_ 000Ig 

003 j 7+0.0005 

84 4 +2 1 

OH- -0.8 


85.06!«;J1 


0.2782 + °™s 

-5 4S99+ 00037 
J ' J -0 0O36 

03 1 3+00001 
87.0+Jo 


km s _1 
km s _1 
AU 

degrees 


Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 


M 

R„ 
logg. 
P. 


i 94+0.06 

J 31+0.05 

4 30+^ 
-Q03 
55 +li " 

"• J -0.05 




1 20 +001 

1 - z ' u -0.01 

j 2i+ 04 
4 33^:8^ 

067^ 

"• u '-0.06 


M 

R G 
[cgs] 

Pg 


Planet mass 

Planet radius 

Planetary surface gravity 

Planet density 

Planet temp 


M„ 

Rp 

Pp 

Teql 


1 76+ 008 

1 -'"-0.14 
i t 1+0.07 
1 - J1 -Q.U 
t ^7+0.09 

• -0.04 
0.78!g 

I960!" 


1 76+ 008 

i 9Q+0.05 
1 - z - 7 -0.12 

3 49+0.O6 

-0.04 

0.82!« : » 


1 90 +0 10 

1 -' u -0.09 

1 20+ 005 
3 47+0.03 

J -^'-0.04 

1.08i| 
1920+j 2 


M, 
R/ 
[cgs] 

Pj 
K 



Table 2. Values for projected stellar rotation rate, spin-orbit misalignment angle and width of the intrinsic stellar line profile com- 
pared to those from the two pr evious Rossiter-Mc Laughlin studies of WASP-3b. Note: The value obtained for v sin i through spectral 
analysis is 13.4 + 1.5 km s" 1 dPollacco et alJl2008l) . 



Parameter Simpson et al. (2009) Tripathi et al. (2010) This study Units 



vsin;' 15.T_\ A 3 I4.l + _\l 13.9^°| km s" 1 

A 13^' 3.3+2-| 5 +6" degrees 



sis we obtained a value for the stellar effective temperature of 
T e ff = 6332 + 105 K. Com paring this temperat ure range with 
the evolutionary tracks from Girardi et al. (2000j) we can get a 
separate mass estimate for the star of M, = 1 .23 ± 0.04 M . This 
is in good agreement with our previous value of M* = 1 .24^q ( ^ 
M . From the maximum error range we can put an upper limit 
on the stellar age of around 2 Gyr. This is an improvement on 
the value presented in lPollacco et al.l J2008) which suggested an 
upper limit of 3.5 Gyr. We cannot put a lower limit on the age of 
the star using Figure 3 as the lower range of the errors lie on or 
beyond the Zero Age Main Sequence. The evolutionary traks of 
Siess et ail (12000) suggest a pre-main sequence lifetime of 0.13 
Gyr for a 1 .2 M star. WASP-3 is a main sequence star, so we can 
impose this value as an extreme lower age limit. In the discovery 
paper a lower age estimate of 0.7 Gyr is presented. We can also 
determine an age estimate from gyrochronology. If we assume 
that the spin axis of WASP-3 is nearly perpendicular to the line 
of sight, then v sin i = 13.9 km s~' and R s = 1.21/? yields a spin 
period of 4.5 sin i days. The period will be shorter than 4.5 days 
if the inclination is substantially less than 90 degrees. 



If we assume that WASP-3 is spinning down because of an- 
gular momentum loss in a hot, magnetically-channeled stellar 
wind, its spin period should increase with time as the square root 
of its age. From their calibration of the period-colour relation in 
Hyades and Praesepe, Delorme et al (2010) find 



625 



11.401 + 12.652(7- 0.631) 



Myr, 



yielding a gyrochronological age of 300 Myr for WASP-3, 
with an uncertainty of order 10 percent. We caution, however, 
that this relation is only applicable once the spin rate has con- 
verged to the asymptotic period-colour relation. While conver- 
gence is probably complete for most stars of this mass by the 
age of 300 Myr, there remains a small possibility that the star 
could have been born as a relatively slow rotator, leading to over- 
estimation of the gyrochronological age. We can, however, state 
with confidence that the gyrochronological age of WASP-3 is 
substantially less than that of the Hyades. 



WASP-3 has 2MASS colour J-K = 0.242. The fastest ro- 
tators at the same colour in the 590-M yr-old Coma Berenices 
open c luster have periods of order 6 days dCollier Cameron et alj 
2009b). Hyades stars of similar c olour also have pe riods in the 
range 5-6 days at age 625 Myr dRadick et al J 1 19871) . as is also 
found to be the case in the Praesepe cluster (age 580 Myr) by 
Delorme et al (2010, MNRAS, submitted). 



4. Conclusions 

The Rossiter-McLaughlin effect present in the WASP-3 sys- 
tem was analysed from observations made using the SOPHIE 
spectrograph on the 1.93m telescope at the Observatoire de 
Haute-Provence (Simpson et al. 2010). We analysed the obser- 
vations using a new method developed bv lCollier Cameron et al.l 
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Fig. 3. The new position of WASP-3 in the T e /f vs R/M 1 ^ 1 plane. The larger of the two boxes shows the data from the discovery 
paper. The smaller box shows the error range from this study. The lines show evolutionary tracks from Girardi et al. (2000) for 1 .2, 
1.3 and 1.4 solar masses and isochrones at 0.1, 0.5, 1, 1.5, 2, 2.5 and 3 billion years. The tracks and isochrones here are for stars 
with solar metalicity, [M/H] = 0. 



(2009a) which involves decomposing the CCF into its various 
components and directly tracking the trajectory of the missing 
starlight across the line profile. This method was incorporated 
into an MCMC analysis of all photometric and spectroscopic 
data available for the WASP-3 system and was found to produce 
good results for the system parameters, helping further constrain 
the values of the spin-orbit misalignment angle, projected stellar 
rotation rate and the impact parameter. The value we obtained for 
the projected spin-orbit misalignment angle A = 9 + c° is close to 
zero a nd agree s with the p r evious results found by Simpso n et al.l 
d2010h and iTripathi et all (l2010l) . Our value of v sin i = 13.9 + 
0.03 km s is in agreement with the value obtained from the 
spectroscopic broadening but determined to a much higher level 
of precision. We conclude that this new method of analysing the 
Rossiter-McLaughlin effect successfully retrieves a more accu- 
rate and precise value for the projected stellar spin rate compared 
with previous measurements, and in doing so it is not vulnerable 
to the systematic error present in previous methods that require 
fitting of Gaussians to non-Gaussian CCFs in order to calculate 
the velocity shifts. It also finds a value for the spin-orbit mis- 
alignment angle that agrees with all previous measurements and 
is of a similar precision. The fact that we clearly detect the signa- 
ture of the missing light after subtracting the model out-of-transit 
profile shows that the method works well for a stellar rotation 
rate that is not much greater th an the intrinsic line width. In fact, 
ICollier Cameron et al.1 (l2009al) showed that this method can be 
successful when v sin i is as low as half the value of the intrinsic 
stellar line width. 



The value we obtained for the spin-orbit alignment angle 
is small enough to support the theory of migration through 
tidal interaction with a protoplanetary dis c. However, recent dis- 
coveries of highly misaligned systems dAnderson et al.l 12009; 
Narit a et al. 2009) suggest that migration cannot be explained 
in all cases purely b y tidal interactions in a protoplanetary disc. 
iHebrard et al.l d2010l) showed that planets with measured spin- 
orbit angles could be sorted into three distinct populations: 
firstly, the majority of hot jupiters that are aligned; secondly, the 
few that are strongly misaligned; and thirdly, the massive plan- 
ets that are mostly moderately but significantly misaligned. They 
propose that this could be the signature of three distinct evolu- 
tion scenarios. Measuring the distribution of alignment angles in 
Hot Jupiter systems and will help inform theories of planetary 
migration. 



Many transit candidates are rejected from RV follow-up 
observations if they are found to be rotating too fast for 
high-precision determination of the RV shift. This is espe- 
cially true for stars of spectral type earlier than F5 which 
tend to be line-poor rapid ro tators. It has already been shown 
(ICollier Cameron et al.l 12010) that using this new method we 
can successfully confirm the existence of a planet-sized object 
by measuring the size of the 'bump' on the stellar line profile, 
opening up a way for establishing the existence of close-orbiting 
planets around rapidly-rotating stars previously inaccessible to 
planet hunters using the RV method. With this new information 
we will also be able to investigate any differences in planet for- 
mation between solar-type and early-type stars. 
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Fig. 4. Mass vs Radius plot for known transiting planets. The two black dots indicate the previous position of WASP-3b from the 
results in the discovery paper and its updated position from the mass and radius values determined in this study. 
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Fig. 2. Top: Residual map of time series CCFs with the model 
spectrum subtracted leaving the bright time-variable feature due 
the light blocked by the planet. Bottom: Here the best-fit model 
for the time-variable feature has also been removed to show the 
overall residual. The horizontal dotted line marks the phase of 
mid-transit. The shift of this line from zero shows the underly- 
ing systemic radial velocity. The two vertical dashed lines are at 
±v sin i from the stellar radial velocity (marked by the vertical 
dotted line). The crosses on the vertical dotted line denote the 
two points of contact at both ingress and egress. 
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